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Abstract 

*The  problem  of  detubling  a  freely  spinning  and  precessing  axi- 
symmetric  satellite  is  considered.  Detumbling  is  achieved  with  another 
axisymmetric  orbital  maneuvering  vehicle  (OMV)  joined  to  the  target  sat¬ 
ellite  with  a  universal  joint.  The  joint  provides  two  rotational  degrees 
of  freedom  and  is  translated  across  the  surface  of  the  OMV  during  the  de¬ 
tumbling  process.  The  target  satellite  and  the  OMV  with  its  three  mo¬ 
mentum  wheels  are  modelled  as  a  five  body  system  using  Eulerian-based 
equations  of  motion  developed  by  Hooker  and  Margulies.  A  Liapunov  tech¬ 
nique  is  applied  to  derive  a  nonlinear  feedback  control  law  which  drives 
the  system  asymtotically  to  a  final  spin-stabilized  state.  State  and 
control  histories  are  presented  and  indicate  that  the  dtumbling  process 
is  benign.  Constraint  force  and  moment  loads  at  the  connection  between 
the  OMV  and  target  satellites  are  also  presented,  and  indicate  that  no 
extreme  loads  are  encountered  during  the  despinning  and  detumbling  process 


Notation 


mass  of  body  A 
total  system  mass 

inertia  dyadic  of  body  A  about  its  mass  center 
angular  velocity  of  body  A 
non-gravitational  external  force  on  body  A 
non-gravitational  external  torque  on  body  A 
geocentric  position  vector  for  mass  center  of  body  A 
interaction  force  acting  on  body  A  through  joint  j 
constraint  torque  acting  on  body  A  through  joint  j 
unit  dyadic 

gravitational  constant 

geocentric  position  vector  for  system  center  of  mass 

unit  vector  in  direction  of  p 

spring-damper  torque  acting  on  body  A  at  joint  j 

unit  vector  along  rotation  axis  of  joint  i 

angle  of  rotation  about  axis  gi 

angular  velocity  of  the  reference  body  (the  OMV) 


THE  DETUMBLING  OF  AN  AXIALLY  SYMMETRIC  SATELLITE  WITH  AN 
ORBITAL  MANEUVERING  VEHICLE  BY  NONLINEAR  FEEDBACK  CONTROL 

I.  Introduction 

The  service  or  repair  of  orbiting  satellites  beyond  direct  reach 
of  the  Space  Shuttle  may  require  an  orbital  maneuvering  vehicle  (OMV)  to 
rendezvous  and  dock  with  the  target  satellites.  If  the  target  -s  spin- 
stabilized  it  may  be  necessary  to  despin  it.  If  the  target  has  experi¬ 
enced  control  system  malfunction  or  for  some  reason  is  not  in  pure  spin, 
it  may  be  ecessary  to  detumble  it.  Docking  followed  by  despinning  or 
detumbling  is  defined  here  as  capture.  Docking  is  accomplished  bv  first 
driving  a  grappling  device  on  the  OMV  to  a  state  of  rest  relative  to  some 
docking  point  on  the  target.  The  OMV  and  target  can  the"  be  connected. 
Despinning  or  detumbling  is  accomplished  by  applying  torques  to  the  tar¬ 
get  through  the  connecting  joint  while  firing  the  OMV  thrusters  to  control 
the  absolute  motion  of  the  two-body  system.  Widhalm  and  Conway  derived  a 
feedback  control  law  (1)  which  solved  the  despinning/detumbling  problem 
for  the  case  of  axisymmetric  target  and  OMV  satellites.  They  used  a  con¬ 
necting  joint  which  could  translate  across  the  surface  of  the  OMV.  The 
translational  degree  of  freedom  of  the  joint  is  depicted  in  Fig.  1  by  the 
double  arrow.  The  ability  to  translate  the  joint  provides  for  joint  pos¬ 
ition  adjustment  during  docking,  and  allows  the  joint  to  be  driven  to  the 
OMV  axis  of  symmetry  during  detumbling.  The  resulting  configuration  can 
then  be  spin-stabilized. 

This  thesis  extends  the  Widhalm  and  Conway  model  to  include  three 
orthogonal  momentum  wheels  on  the  OMV,  and  develops  a  feedback  control 


law  to  couple  the  momentum  wheel  torques  to  the  system  state.  The  target 
satellite  docked  with  the  OMV  and  its  three  momentum  wheels  are  modelled 
as  a  system  of  five  rigid,  constant  mass  bodies.  The  control  problem  is 
formulated  by  defining  the  system  initial  configuration  and  the  desired 
final  state,  and  by  deriving  the  equations  of  motion. 


II.  Problem  Formulation 


System  Configuration 

An  axisymmetric  target  satellite  is  docked  with  an  OMV  which, 
with  its  three  orthogonal  momentum  wheels,  is  also  axisymmetric  (see 
Fig.  2).  The  target  and  OMV  are  connected  with  a  universal  joint  having 
two  rotational  degrees  of  freedom  and  the  capability  of  translation 
across  the  surface  of  the  OMV.  The  translational  degree  of  freedom  is 
depicted  in  Fig.  2  by  the  bold  double  arrow.  The  center  of  mass  of  the 
target  lies  on  the  83  axis  as  does  the  mass  center  of  the  OMV-momentum 
wheel  combination.  The  OMV  is  in  a  state  of  pure  spin  about  63,  and  the 
target  is  in  a  state  of  spin  with  precession  about  63  at  a  rate  equal  to 
the  OMV  spin  rate.  With  no  external  moments  or  forces  acting  on  the  sys¬ 
tem  a  dynamically  stable  configuration  results.  This  configuration  rep¬ 
resents  the  initial  state  of  the  system.  The  detumbling  and  despinning 
process  is  complete  when  the  joint  has  been  driven  to  a  position  lying  on 
the  63  axis,  and  the  target  spin  rate  relative  to  the  OMV  is  zero  with 
the  OMV  itself  still  in  a  state  of  pure  spin  about  the  83  axis.  This 
configuration,  or  one  arbitrarily  close  to  it,  represents  the  desired 
final  state  of  the  system.  Both  the  initial  and  final  states  as  defined 
imply  that  the  initial  and  final  angular  velocities  of  the  momentum 
wheels  having  rotational  freedom  about  the  8^  and  bg  axes  are  zero.  The 
initial  and  final  angular  velocity  of  the  83  wheel  is  arbitrary. 

The  Hooker-Margulies  Equations 

The  dynamical  attitude  equations  for  a  two-body  satellite  were 
derived  by  Fletcher,  Rongved  and  Yu  (2),  and  were  generalized  for  an 


complete  derivation  is  a  set  of  3n  scalar  equations  for  an  n-body  system. 
The  equations  are  free  of  the  unknown  joint  constraint  forces,  but  still 
contain  the  unknown  constraint  torques  (3:125). 

Hooker  (4)  showed  in  a  subsequent  paper  that  the  constraint 
torques  could  also  be  explicitly  eliminated.  The  equations  derived  in 
reference  (3)  are  written  for  all  the  bodies  lying  on  one  side  of  a  se¬ 
lected  joint  and  subsequently  added.  The  interaction  torques  all  cancel 
in  pairs,  with  the  exception  of  the  constraint  torque  at  the  free  joint. 

If  the  selected  joint  has  a  rotational  degree  of  freedom  about  an  axis 
g,  the  dot  product  of  g  and  the  expression  just  found  for  the  torque  is 
zero.  Writing  the  dot  product  and  setting  it  to  zero  yields  an  equa¬ 
tion  free  of  the  constraint  torque.  Repeating  the  process  for  each  de¬ 
gree  of  freedom  at  each  joint  eliminates  all  the  unknown  constraint 
torques  and  yields  a  system  of  r  equations  for  an  n-body  system  having 
r  rotational  degrees  of  freedom.  These  equations  are  referred  to  in 
some  of  the  references  as  the  modified  Hooker-Margulies  equations. 

Although  for  an  identical  dynamical  system  a  Lagrangian  deriva¬ 
tion  would  provide  the  same  number  of  equations  as  the  HM  equations,  the 
resulting  expressions  would  not  be  written  in  terms  of  the  physical  body 
axes,  as  are  the  HM  and  modified  HM  equations.  As  a  result  adaptation 
to  active  control  and  modification  to  include  effects  such  as  joint  motion 
would  be  more  difficult  (4:1205).  The  modified  HM  equations  significant¬ 
ly  reduce  the  required  computer  time  for  solution  due  to  the  reduction 
in  the  number  of  equations  from  3n  to  r. 

The  constraint  torques  at  the  joints  were  eliminated  explicitly 
by  taking  the  dot  product  of  the  unit  vector  about  which  the  joint  is 


n-body  system  by  Hooker  and  Margulies  (3).  Both  derivations  assume  that 
the  bodies  are  connected  by  joints  which  are  fixed  with  respect  to  the 
bodies  they  connect.  In  addition,  the  generalized  Hooker-Margulies  ( HM) 
equations  assume  that  the  bodies  are  '■onnected  in  a  topological  tree. 

This  means  that  there  are  no  closed  loops  formed  by  the  interconnected 
bodies.  The  restriction  of  immoveable  joint  is  removed  in  an  extension 
of  the  HM  equations  which  will  be  discussed  later.  The  extended  equa¬ 
tions  then  apply  to  the  OMV-target  satellite  system  when  modelled  as  a 
system  of  five  interconnected  rigid  bodies. 

The  derivation  of  the  HM  equations  begins  with  the  Newton  and 
Euler  equations  for  an  n-body  system.  Each  of  the  two  sets  of  equations 
contains  force  terms  representing  the  unknown  constraint  reactions  which 
occur  at  the  joints  between  adjacent  bodies.  An  expression  for  each 
joint  constraint  force  can  be  isolated  by  writing  Newton's  equations  for 
all  the  bodies  that  lie  to  one  side  of  any  selected  joint.  The  equations 
are  added  together,  with  the  result  that  all  the  interaction  forces  can- 
in  pairs  with  the  exception  of  the  constraint  force  occurring  at  the 
selected  joint.  Repeating  the  process  for  all  the  joints  in  the  system 
yields  expressions  for  all  the  unknown  constraints  in  terms  of  the  system 
external  forces  and  in  terms  of  the  inertial  accelerations  of  each  body 
in  the  system. 

The  joint  constraint  forces  appear  in  Euler's  equations  as 
torques  about  the  individual  body  mass  centers,  and  can  be  replaced  by 
the  expressions  for  the  constraints  obtained  in  the  process  described. 

The  result  is  the  original  Euler  equations  for  the  system,  but  with  the 
unknown  joint  interactions  explicitly  eliminated.  The  result  of  the 


free  to  rotate  (there  may  of  course  be  more  than  just  one)  with  the 
summed  vector  dynamical  equations  for  the  bodies  that  lie  to  one  side  of 
that  joint.  The  value  of  the  constraint  torque  can  be  computed  after 
the  modified  HM  equations  are  solved,  however.  Given  the  system  state, 
the  derivatives  of  the  state  variables  can  be  computed  using  the  modi¬ 
fied  HM  equations.  All  the  variables  in  the  equation  which  was  dotted 
with  the  joint  degree  of  freedom  vector  g  are  then  known,  and  can  all  be 
brought  to  one  side  of  the  equation  with  the  resulting  sum  equalling  the 
scalar  components  of  the  constraint  torque  (4:1207). 

The  interaction  force  at  the  joints  can  be  computed  by  finding 
the  accleration  of  the  mass  center  of  the  system  of  bodies  lying  to  one 
side  of  a  joint  (relative  to  the  system  mass  center)  and  multiplying  by 
the  total  mass  of  that  subsystem.  For  the  case  of  no  external  forces, 
the  product  is  the  vector  constraint  force  acting  at  the  selected  joint. 


Application  of  the  HM  Equations 

The  Eulerian-based  equations  of  motion  for  multi-body  systems 
given  by  Hooker  and  Margulies  (3)  and  modified  by  Hooker  (4)  are  restric¬ 
ted  to  those  systems  of  bodies  connected  in  such  a  way  that  no  closed 
loops  are  formed,  and  do  not  account  for  motion  of  the  joints  relative 
to  the  bodies  adjacent  to  the  joint.  This  last  restriction  was  removed 
in  an  extenstion  of  the  modified  HM  equations  by  Conway  and  Widhalm  (5)  to 
permit  the  translation  of  the  joint  across  the  face  of  the  OMV.  Thus  the 
extended  equations  can  be  applied  to  the  OMV,  momentum  wheels  and  target 
satellite  system  under  consideration.  Referring  to  Fig.  2  the  OMV  is 
labelled  body  0,  with  the  b  basis  fixed  at  its  geometric  center.  The 
geometric  center  of  the  OMV  is  chosen  to  coincide  with  the  center  of  mass 


of  the  OMV -momentum  wheel  composite.  That  is,  remove  the  target  from  the 

system  and  the  system  center  of  mass  then  lies  at  the  origin  of  the  b 

frame.  The  momentum  wheels  lying  on  axes  6^,  b^,  and  b^  are  labelled 

bodies  1,  2,  and  3  respectively.  The  target  satellite  is  body  4.  The 

e  basis  is  fixed  in  body  0  at  its  mass  center,  and  the  n  basis  is  fixed 

in  the  target  at  its  mass  center.  Since  the  target  is  precessing  about 

its  own  angular  momentum  vector  at  a  rate,  ip  ,  the  OMV  is  positioned 

relative  to  the  target  so  that  the  target's  center  of  mass  and  angular 

momentum  vector  both  lie  on  the  b^  axis.  The  OMV  is  then  spun  about  bj 

at  the  same  rate,  ip  ,  the  target  precession  rate.  The  cone  angle, 

and  the  distance  from  the  target  center  of  mass  to  the  joint  determine 

the  required  position  of  the  joint  on  the  OMV  face  for  docking.  The 

•  • 

target's  cone  angle,  ,  precession  rate,  \p  ,  spin  rate,X5,  and  mass 
properties  are  related  as  shown  by  Greenwood  (6:386)  and  repeated  here: 


P  =  IV{I0  '  I)cosy4 


(1) 


where  I  is  the  target's  moment  of  inertia  about  the  n^  and  axes  and 
IQ  is  the  target's  moment  of  inertia  about  the  spin  axis,  n^. 

The  two  rotational  degrees  of  freedom  required  at  the  universal 
joint  include  rotation,^,  about  an  axis  g4  parallel  to  b1  and  rotation 
y5  ,  about  an  axis  g5  parallel  to  n^.  The  rotational  degrees  of  free¬ 
dom  for  the  momentum  wheels  labelled  1,  2,  and  3  are  axes  g^,  gj,  and  g^ 
respectively.  These  three  axes  are  parallel  to  the  corresponding  6  frame 
axes,  with  the  wheel  joints  themselves  taken  to  be  at  the  mass  centers  of 
the  wheels.  From  this  initial  docked  configuration  the  problem  is  to 
drive  the  system  to  a  final  spin-stabilized  state  with  a  set  of  feedba 
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controls.  This  final  state  is  specified  by  requiring  the  joint  location 
to  coincide  with  the  axis,  and  that  the  cone  angle,  ,  and  spin  rate 
y?  ,  be  reduced  to  steady  state  values  of  zero. 

In  the  following  equations  of  motion,  all  vectors  and  scalar 
rates  are  with  respect  to  the  e  basis  fixed  in  the  main  body,  the  OMV. 

The  attitude  equations  for  the  OMV-target  satellite  system  can  be  derived 
directly  from  the  extended  equations  given  by  Conway  and  Widhalm  (5)  and 
are: 
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1,  if  gi  belongs  to  a  joint  anywhere  on  the  chain 
of  bodies  connecting  n  and  the  reference  body  (0) 
0,  otherwise  (e.g.  if  0) 
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The  vector  is  the  vector  from  the  center  of  mass  of  body  A 
to  the  joint  leading  to  body  fi.  E^  is  determined  from 

i‘*„Kk  <91 

and  Ex  is  the  vector 

h  -  3^'3  *  *  *\x-  *  ■  SA*  *AA'  «>A  ♦  T;  +  S  T9° 

J  ^X 

+  h  x  fa+„|xVx  cl>+m£v  (V  W  no) 

+  m p"3(  1  -  3?0)  .  DjjJ 

In  Eq  (2)  <DQ  =  Wq-^  +  <*\)2®2  +w03®3»  and  the  ^i  correspond 
to  degrees  of  freedom  about  the  unit  vectors  g^ .  Superscript  R  implies 
that  the  indicated  time  differentiation  is  performed  with  respect  to  an 
observer  fixed  in  the  e-frame.  From  Eqs  (7)  and  (8)  and  from  the  defi¬ 
nition  given  for  L it  is  clear  that  the  indicated  time  derivatives 
of  Dq4  are  determined  from  the  universal  joint  motion  alone  since  the 
momentum  wheel  joints  do  not  move  with  respect  to  the  OMV.  The  joint 
motion  is  specified  as  are  the  control  torques  T3^  in  Eq  (10).  In  this 
analysis,  the  external  torques  T'  are  the  three  orthogonal  thruster 
torques  acting  on  the  OMV,  and  all  external  forces  are  assumed  zero.  In 
addition,  all  gravitational  terms  are  ignored  (all  terms  containing  the 
position  vector  p  ).  The  way  in  which  the  control  torques  (including 


& 

thruster  torques)  and  the  joint  motion  is  specified  is  covered  in  the 
next  chapter. 

The  constraint  forces  and  torques  have  been  eliminated  from  the 
equations  of  motion  as  described  earlier.  These  quantities  can  be  deter¬ 
mined  by  the  methods  described  in  references  (3)  and  (4).  The  expres¬ 
sions  for  the  constraint  force  and  torque  at  the  universal  joint  are: 

F0J  =  (m  -  m4)  fQ  (11) 

Ln  =  (  *44  +  04O5  ‘  ^0  +  (  *44  *  V  ^4 

..  _  (12) 

+  ^  ^44  •  9g  )  y  5  ^4  m®40  x  ® 

where  ?Q  is  the  position  vector  of  the  mass  center  of  the  OMV-momentum 
wheel  combination,  relative  to  the  system  center  of  mass,  and  where  G  is 


III.  Nonlinear  Feedback  Control 


The  detumbling  and  despinning  problem  presented  here  involves 
driving  the  universal  joint  connecting  the  OMV  and  target  satellites 
across  the  face  of  the  OMV  to  a  point  coincident  with  the  bj  axis  (see 
Fig.  2).  Feedback  control  is  used  to  maintain  the  attitude  of  the  five- 
body  system  in  such  a  way  that  the  final  state  of  the  system  is  spin 
stabilized.  The  Droblem  solution  incorporates  a  feedback  control  ap¬ 
proach  in  which  an  eight  element  control  vector,  u,  is  a  nonlinear  func¬ 
tion  of  the  system  state  variables.  In  this  chapter  the  eight  controls 
are  defined  so  that  the  equations  of  motion  can  be  written  in  a  simpler 
form  more  suitable  for  Liapunov  analysis.  This  procedure  closely  paral¬ 
lels  that  described  by  Widhalm  and  Conway  (7:6-9)  in  their  derivation  of 
a  control  law  for  the  two-body  satellite  system  described  earlier. 
Liapunov’s  direct  method  is  then  used  to  derive  a  control  law  which  is 
globally  asymptotically  stable  with  respect  to  the  spin-stabilized  equi¬ 
librium  state. 

Eq  (2)  is  first  written  in  the  form 

A  £  =  F*  (13 

where  A  is  defined  as  the  8  x  8  matrix  on  the  left-hand  side  of  Eq  (2). 
The  vector,  F*.  is  defined  as  the  eight  element  vector  on  the  right-hand 


Eq  (14)  is  derived  from  the  state  variables 


x  =  [  x]  x2  x3  x4  x5  xg  x7  x8] 


=  [  W01  w02  w03  *1  *2  *3  *4  *5  ^ 


(15) 


The  control  vector,  u,  can  be  selected  (7:6)  so  that  Eq  (13)  can 
be  written  as 


A  7  =  F  +  u  (16) 

Since  the  components  of  the  vector,  7,  are  the  three  angular  acceleration 
components  of  the  OMV,  the  three  scalar  angular  accelerations  of  the 
othogonal  momentum  wheels,  and  the  angular  accelerations  of  the  target 
about  the  two  degrees  of  freedom  at  the  joint,  the  appropriate  control 
vector  components  are  apparent.  The  control  vector,  u,  is  composed  of 
three  orthogonal  (thruster)  torques  about  each  of  the  e-basis  vectors, 
three  internal  torques  applied  at  the  wheel  axes,  and  finally  two  inter¬ 
nal  torques  applied  in  the  two  degrees  of  freedom  of  the  universal  joint. 
These  control  torques  are  designated  u^  through  ug,  respectively.  Pre¬ 
multiplying  Eq  (16)  by  the  inverse  of  matrix  A  yields 

7  =  A-1  F  +  A-1  u  (17) 

since  the  matrix.  A,  will  always  be  invertible  for  physical  systems.  The 
system  of  Eq  (2)  is  augmented  with  the  kinematical  equation 

Xg  =  Xj  (18) 

where  Xg  is  defined  as  the  target  precession  angle,  7^,  and  competes 
the  set  of  equations  of  attitude  motion.  We  define  the  augmented  state 


vector,  .jXg,  to  contain  the  vector,  "x,  plus  the  ninth  element,  Xg. 

To  derive  the  control  vector,  u,  as  a  nonlinear  function  of  the 
augmented  state  vector,  ^Xg,  and  the  joint  motion,  a  lemma  presented  by 
Vidyasagar  (8)  is  applied.  The  lemma  applies  to  autonomous  systems,  and 
the  system  of  Eq  (16)  is  nonautonomous  because  the  OMV-target  connecting 
joint  motion  is  a  specified  function  of  time.  Widhalm  and  Conway  (7:7) 
have  suggested  that,  by  specifying  the  joint  motion  as  a  third  order  lin¬ 
ear  system  which  is  asymptotically  stable  with  respect  to  the  desired 
final  joint  position,  the  third  order  system  with  Eqs  (17)  and  (18)  form 
an  autonomous  system.  Since  the  desired  final  joint  position  relative 
to  the  B-frame  is  given  by  the  vector  [0  0  c  ]  ,  and  since  the  HM 
equations  are  written  in  terms  of  the  e-frame,  we  define  the  vector,  c, 
to  be  the  vector  leading  from  the  origin  of  the  b  basis  to  the  mass 
center  of  the  OMV,  the  origin  of  the  e-frame.  The  scalar  values  that  are 
to  be  driven  to  zero  for  asymptotically  stable  joint  motion  are  the  posi¬ 
tion,  velocity,  and  acceleration  components  of  the  joint  lying  along  the 
62  axis.  These  scalars  are  given  by: 

(  c  +  Lo4  )  •  b2  =  y] 

*-04  **2  =  y2  =  ^1 

^04  '  *2  =  y3  =  h 

The  derivative  of  the  vector,  c,  does  not  appear  since  the  vector  posi¬ 
tion  of  the  OMV  mass  center  relative  to  the  geometric  center  of  the  OMV 
is  constant.  The  complete  autonomous  system  required  for  the  application 
of  the  previously  mentioned  lemma  is  given  by  Eqs  (17),  (18),  and  (19) 


and  is  formed  by  writing 


•  -l  -  -1 

x  =  A  1  F  +  A  1  u 

Xg  =  x7  (20) 

j  =  Dy 

where  D  is  a  negative  definite  matrix  selected  to  obtain  the  desired  de¬ 
cay  of  the  scalar  y^  given  in  Eq  (19).  The  joint  motion  terms  contained 
in  the  matrix.  A,  and  in  the  vector,  F,  are  now  specified  by  the  vector 
y  leaving  the  system,  (20),  independent  of  time. 

The  lemma  developed  by  Vidyasagar  (8:157)  in  his  discussion  of 
Liapunov's  direct  method  now  applies  to  the  system,  (20),  and  is  stated 
as  follows:  "Let  V(^Xg,  y)  be  continuously  differentiable  and  suppose 
that  for  some  d  >  0  the  set 

=  [  ^g,  y  :  V ( i"Xg ,  y)  <  d  ] 

is  bounded.  Suppose  that  V  is  bounded  below  over  the  set  and  that 
V(jXg,  y)  ^  0  for  all  jXg  and  y  in  S^.  Let  S  denote  the  subset  of  de¬ 
fined  by 

S  =  [  iXg,  y  :  V(-|Xg,  y)  =  0  ] 

and  let  M  be  the  largest  invariant  set  of  a  system  which  is  contained  in 
S.  Then  whenever  ^Xg(O)  and  y(0)  are  members  of  S^,  the  solution  of  the 
system,  (20),  approaches  M  as  t  — 

Now  M  is  an  invariant  set  of  system,  (20),  if  every  trajectory 
starting  from  an  initial  point  in  M  remains  in  M  for  all  time.  Since  the 
system,  (20),  is  autonomous  every  trajectory  through  its  state  space  is 
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an  invariant  set  (8:156).  The  task  then  is  to  find  a  candidate  Liapunov 
function,  V,  that  meets  the  requirements  of  the  lemma.  To  derive  a  non¬ 
linear  feedback  control  law  that  drives  the  five-body  system  to  the  de¬ 
sired  final  state  of  spin-stabilized  equilibrium,  Widhalm  and  Conway  (7) 
suggested  a  candidate  Liapunov  function  like 

V  =  (1/2)  xT  I  x  +  (1/2)  K  x*  +  yT  R  y  (21) 

where  I  is  the  identity  matrix,  K  is  a  positive  constant,  and  R  is  a 
positive  definite  constant  matrix.  The  function  is  continuously  differ¬ 
entiable,  and  it  is  easy  to  select  a  vector  ^xg  and  a  constant,  d,  to 
demonstrate  that  a  set  like  S*j  exists. 

The  condition  on  V(jxg,  y)  must  be  satisfied;  differentiating 
V  with  respect  to  time  yields 

V  =  xT  I  7.  +  K  Xg  xg  +  R  y  +  yT  R  y  (22) 

Substitution  from  the  system,  (20),  gives  the  result 


V  =  xT  I  [  A'1  F  +  A'1  u  ]  +  K  x9  x7 

+  yT  [  DT  R  +  R  D  ]  y  (23) 

Now  since  R  is  a  positive  definite  matrix  and  D  was  specified  to  be  a 
negative  definite  matrix*  then  the  expression  DTR  +  RD  is  negative 
definite.  Writing  the  Liapunov  matrix  equation 

-Q  =  0T  R  +  R  D  (24) 

and  defining 

V*  =  -  yT  Q  yT  (25) 
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implies  both  that  Q  is  positive  definite  from  Eq  (24)  and  that  as  a 
★  •  ★ 

result  both  V  and  -V  are  positive  definite  (8:172).  Thus  the  third 
term  in  Eq  (23)  is  negative  definite.  To  make  V  at  least  negative  semi- 
definite,  select  the  control  vector 

u  =  -F  +  A  [  0  0  0  0  0  0  -Kx9  0  ]T  -  A  B  x  (26) 

where  the  matrix  B  is  positive  definite  or  positive  semi-definite.  If 
the  matrix  B  is  positive  definite,  then  V  is  negative  semi -definite  in 
Xg.  However,  if  the  elements  of  B  are  positive  except  for  Bgg  =  0  , 

V  is  negative  semi-definite  in  Xg  and  Xg.  Then  from  the  above  lemma  the 
system,  (20),  with  control  vector  u  of  Eq  (26)  is  asymptotically  stable 
with  respect  to  the  largest  invariant  set  contained  in  the  Xg,  Xg  plane. 
But  from  Eq  (2b)  it  can  be  seen  that  any  non-zero  xg  results  in  a  non¬ 
zero  control  vector  u,  which  in  turn  causes  a  departure  from  the  Xg,  xg 
plane.  Therefore  the  largest  invariant  set  lies  in  the  Xg,  xg  plane,  but 
only  in  that  region  where  Xg  is  zero,  namely  the  Xg  axis.  The  control 
vector  u  of  Eq  (26)  is  then  the  desired  nonlinear  feedback  control  law 
for  spin-stabilization  of  this  system. 

Non- zero  off-diagonal  terms  must  be  included  in  the  B  matrix  of 
Eq  (26)  which  will  couple  the  momentum  wheel  torques  to  state  variables 
other  than  just  the  momentum  wheel  angular  velocities.  The  torques  must 
be  coupled  to  the  state  variables  in  such  a  way  that  the  angular  rates  of 
the  g^  and  g2  wheels  decay  to  zero,  so  that  the  system  stabilizes  in  pure 
spin.  Implementing  the  control  law  of  Eq  (26)  requires  the  determination 
of  the  non-zero  elements  of  the  matrix  B,  and  the  constant  K.  In  their 
work  with  the  two-body  satellite  system,  Widhalm  and  Conway  have  suggested 


that  controls  can  be  kept  within  reasonable  limits  by  ensuring  that  the 
target  center  of  mass  does  not  depart  appreciably  from  the  bg  axis.  This 
can  be  accomplished  by  selection  of  the  matrix  D,  which  controls  the  uni¬ 
versal  joint  motion,  in  conjunction  with  K  and  the  Byy  element  of  the  B 
matrix.  Joint  motion  can  be  specified  which  will  closely  follow  the  de¬ 
cay  of  the  target's  precession  angle,  and  maintain  the  proximity  of  the 
target  mass  center  to  the  bg  spin  axis.  Finally,  substitution  of  the 
control  vector  u  of  Eq  (26)  into  Eq  (20)  yields  the  complete,  linear 
system 

7=[000000  -Kxg  0  ]T  -  B  x 
Xg  =  Xy  (27) 

y  =  D  y 

From  Eq  (27)  it  can  be  clearly  seen  that  if  the  matrix  B  is  dia¬ 
gonal,  and  if  the  initial  momentum  wheel  angular  rates  as  well  as  the 
OMV  angular  velocity  components  and  wQ2  are  zero,  they  will  remain 

zero  throughout  the  detumbling  process.  In  the  case  where  the  only  off- 
diagonal  terms  in  B  are  those  coupling  the  momentum  wheel  torques  to 
selected  system  states,  the  a»01  and  wQ2  angular  velocity  components  of 
the  OMV  still  remain  zero.  This  implies  that  the  throughout  the  maneuver 
the  OMV  remains  in  a  state  of  pure  spin. 


IV.  Results 
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In  this  chapter  the  basis  for  the  selection  of  initial  conditions 
and  system  mass  properties  is  presented,  as  well  as  the  values  selected. 
System  state  histories  and  control  torque  histories  are  presented  for 
four  detumbling  maneuvers.  The  four  maneuvers  include  detumbling  with 
feedback  control  without  momentum  wheel  coupling,  feedback  control  with 
a  single  momentum  wheel  coupled,  detumbling  using  only  gimbal  (joint) 
torques,  and  finally,  detumbling  using  only  gimbal  torques  for  the  first 
250  seconds,  followed  by  50  seconds  of  feedback  control  applying  the  full 
control  torque  vector  with  one  wheel  coupled.  Representative  histories 
for  the  constraint  loads  at  the  joint  between  the  OMV  and  the  target  sat¬ 
ellite  are  also  given. 

System  mass  properties  were  selected  based  on  values  used  in  pre¬ 
vious  research  (9).  The  mass  properties  of  the  composite  body  consisting 
of  the  OMV  and  the  three  attached  momentum  wheels  were  selected  to  dupli¬ 
cate  those  of  the  OMV  model  used  by  Widhalm  (9).  Reasonable  but  arbitrary 
values  of  mass  and  moments  of  inertia  were  then  selected  for  the  three 
identical  wheels.  From  this  information  and  by  specifying  the  location 
of  the  mass  center  of  the  composite  OMV-momentum  wheel  system,  the  mass 
and  inertia  matrix  for  the  OMV  (without  wheels)  were  computed.  In  this 
way  direct  comparisons  could  be  made  between  the  two-body  and  five-body 
analyses  by  completely  decoupling  the  wheel  control  torques  from  all  non¬ 
zero  system  states.  As  mentioned  previously,  this  is  accomplished  by 
selecting  all  off-diagonal  terms  of  the  matrix  B  in  the  linear  system 
(27)  to  be  zero.  Under  this  condition,  the  values  of  elements  B^,  Bgg, 
and  Bgg  are  completely  arbitrary.  In  addition  to  the  various  mass 
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properties,  the  system  state  variable  initial  conditions  were  also  dupli¬ 
cated  from  reference  (9).  The  mass  properties  and  initial  conditions  are 
given  in  Tables  I  and  II,  respectively. 

From  Eq  (27)  and  from  the  discussion  that  immediately  follows  it, 
it  can  be  seen  that  the  values  selected  for  elements  and  B22  are  a^so 
arbitrary,  since  they  will  always  be  multiplied  by  states  and  x^,  which 
remain  equal  to  zero  throughout  the  detumbling,  despinning  maneuver.  The 
values  actually  selected  are  those  used  by  Widhalm  (9)  in  the  application 
of  feedback  control  to  the  two-body  satellite  system.  In  that  case  states 
x1  and  x2  attained  non-zero  values  after  approximately  290  seconds  of 
open-loop  control.  Feedback  control  was  then  applied  in  an  attempt  to 
spin-stabilize  the  system,  the  desired  response  being  obtained  using  the 
values  =  B^  =  0.046  . 

The  remaining  constants  to  be  determined  include  B77>  BQ8,  K,  and 
the  elements  of  the  D  matrix  of  the  third  order  joint  motion  equation  in 
Eq  (27).  Selecting  the  scalar  equation  from  the  12  equation  system  of 
Eq  (27)  corresponding  to  the  target  precession  angle  motion,  the  constants 
B77  and  K  can  be  determined  by  specifying  a  desired  final  precession 
angle  at  the  end  of  the  maneuver.  A  total  maneuver  time  of  300  seconds 
was  selected,  and  critical  damping  specified.  A  final  precession  angle 
corresponding  to  the  target  mass  center  on  the  OMV  spin  axis  and  a  final 
joint  position  equal  to  0.05%  of  initial  joint  position  was  specified. 

The  result  obtained  from  these  requirements  is  a  solution  with  equal 
eigenvalues  of  about  -0.035,  with  B77  =  0.07  and  K  =  0.00123.  In  a 
similar  way,  with  the  specification  of  the  final  joint  position  just 
given  and  with  critical  damping  of  the  joint  motion,  equal  eigenvalues 
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TABLE  I 


Satellite  Mass  Properties 


MASS  Ij  I2 

*3 

Target  Satellite 

1000  Kg  1000  Kg-m2  1000  Kg-m2 

1100  Kg-m2 

OMV 

4500  Kg  6400  Kg-m2  6400  Kg-m2 

11800  Kg-m2 

Momentum  Wheels 

10  Kg  25  Kg-m2  25  Kg-m2 

55  Kg-m2 

TABLE  II 

Initial  Conditions 

Variable 

Meaning 

Value 

W01 

e1  component  OMV  angular  velocity 

0.0 

W02 

@2  component  OMV  angular  velocity 

0.0 

W03 

e.j  component  OMV  angular  velocity 

0.102 

*1 

#1  momentum  wheel  spin  rate 

0.0 

*2 

#2  momentum  wheel  spin  rate 

0.0 

*3 

#3  momentum  wheel  spin  rate 

0.0 

target  precession  angle  rate 

0.0 

*5 

target  spin  rate 

0.009 

(Lq4  +  c)  •  b2 

joint  position 

0.599 

4)4  '  ^2 

joint  velocity 

0.0 

4)4  *  ^2 

joint  acceleration 

0.0 

V4 

target  precession  angle 

0.349 

Note:  Values  in  TABLE  II  are  given  in  meters,  radians, 
or  radians  per  second,  as  applicable. 


of  -0.04  were  obtained  for  the  joint  motion  equation.  This  resulted  in 
values  of  =  -0.000064,  =  -0.0048,  and  Dgg  =  -0.12  .  Finally, 

solving  the  eighth  scalar  equation  of  the  system  of  Eq  (27)  by  requiring 
that  the  target  spin  rate  be  reduced  to  0.5%  of  its  initial  value  yields 
the  value  of  Bgg  =  0.02  . 

With  the  momentum  wheels  uncoupled  the  system  behavior  was  ident¬ 
ical  to  that  of  the  two-body  system  in  Fig.  1,  verified  by  the  comparison 
of  the  results  with  those  of  Widhalm  (9).  Wheel  torques  were  zero,  with 
the  0MV  thrust  torques  and  universal  joint  torques  given  in  Figs.  3  and  4. 
None  of  these  control  torques  reached  large  values,  and  it  can  be  seen 
that  the  controls  vary  smoothly  with  time  with  no  abrupt  changes.  The 
joint  motion  and  precession  angle  decay  behaved  as  specified,  decaying 
to  the  small  final  values  specified  with  no  overshoot.  The  results  are 
given  in  Fig.  5.  The  target  spin  rate  decay  displayed  similar  behavior. 
The  precession  angle  rate  of  change  and  target  spin  rate  histories  are 
given  in  Fig.  6,  where  it  can  be  seen  that  no  radical  motion  has  occurred. 
The  constraint  loads  at  the  universal  joint  due  to  lack  of  rotational 
freedom  and  due  to  the  loading  required  to  drive  the  joint  across  the  0MV 
were  relatively  small.  The  constraint  torque  maximum  value  was  approxi¬ 
mately  2  Nm,  and  the  maximum  force  encountered  in  driving  the  joint  was 
about  1  N.  These  results  are  shown  in  Fig.  7,  where  the  magnitude  of  the 
constraint  force  is  plotted,  as  well  as  the  component  of  that  force  lying 
along  the  axis  of  joint  motion.  This  component  corresponds  to  the  load 
capacity  that  would  be  required  for  a  jackscrew,  for  instance,  were  such 
a  device  used  on  the  0MV  to  obtain  the  desired  joint  motion. 

Although  global  asymptotic  stability  is  only  guaranteed  by  the 


lemma  used  earlier  for  the  case  of  feedback  control  with  the  control  law 
of  Eq  (26),  it  was  of  interest  to  attempt  system  despin  and  detumbling 
using  only  control  torques  applied  at  the  universal  joint.  The  result 
was  that  the  system  displayed  no  radical  or  unusual  behavior,  and  appro¬ 
ached  a  nearly  spin-stabilized  state  after  300  seconds  of  control.  The 
OMV  still  remained  in  a  state  closely  approximating  pure  spin,  although 
non-zero  angular  velocity  components  did  remain.  Figs.  8  and  9  show  the 
OMV  angular  velocity  history  for  this  maneuver,  and  indicate  that  the  de¬ 
parture  from  pure  spin  was  modest.  Constraint  torque  and  interaction 
force  at  the  connection  also  remained  small,  although  as  indicated  in 
Figs.  10  and  11,  a  quasi-steady  state  condition  was  reached  with  continu¬ 
ous  constraint  and  control  torques  experienced  at  the  joint.  A  residual 
precession  angle  remained  after  the  maneuver,  but  the  rate  of  change  of 
the  precession  angle  was  very  small,  indicating  that  the  available  con¬ 
trol  torques  were  able  to  maintain  the  system  configuration  but  unable  to 
change  it  at  any  appreciable  rate.  Although  there  was  some  oscillation 
in  the  precession  angle  rate  of  change,  no  other  violent  behavior  occur- 
ed  with  respect  to  either  target  spin  or  precession  angle  changes.  The 
results  are  shown  in  Figs.  12  and  13. 

The  maneuver  described  above  was  repeated,  but  at  250  seconds  the 
full  control  vector  was  applied  in  an  attempt  to  spin-stabilize  the  sys¬ 
tem.  The  residual  non-zero  c component  of  the  OMV  angular  velocity 
was  reduced  nearly  to  zero,  and  the  OMV  spin  rate  brought  to  a  steady 
value  (see  Figs.  14  and  15).  The  relative  spin  rate  of  the  target  and  the 
precession  angle  rate  of  change  were  likewise  reduced  to  very  nearly  zero 
as  indicated  in  Fig.  16.  Although  not  plotted,  the  precession  angle  of 
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the  target  was  reduced  to  half  of  its  value  at  250  seconds,  or  approxi¬ 
mately  0.9  degrees.  Figs.  17  through  19  show  the  attendant  reductions 
in  the  constraint  and  control  torques  and  in  the  constraint  force  at  the 
joint.  These  data  clearly  indicate  that  the  system  is  moving  toward  a 
spin  stabilized  equilibrium. 

An  attempt  was  made  to  couple  the  momentum  wheels  to  the  target 
precession  angle  rate  and  spin  rate.  The  magnitude  of  the  OMV  thruster 

A 

torque  about  each  of  the  b  axes  was  integrated  over  the  maneuver  period, 
and  the  three  resulting  values  added  to  obtain  some  measure  of  the  total 
torque  power  required  from  the  thrusters.  The  resulting  total  was  then 
compared  for  maneuvers  completed  using  various  off-diagonal  gain  terms  in 

A 

the  B  matrix  of  Eq  (26).  No  attempt  at  using  control  torque  on  the  b^ 
wheel  yielded  a  reduction  in  the  OMV  thruster  torque  requirements.  Using 

A 

applied  torque  at  the  b2  wheel  by  setting  B57  to  1.0  resulted  in  a  reduc¬ 
tion  of  the  integrated  thrust  torque  from  436  Nm-sec  for  no  wheel  coupling 
to  approximately  377  Nm-sec.  The  resulting  thruster  control  torques  and 
momentum  wheel  control  torque  are  shown  in  Figs.  20  and  21.  Fig.  21  also 

A 

includes  the  b2  momentum  wheel  angular  rate  history  during  the  maneuver. 

These  results  indicate  that  the  system  is  well  behaved  in  the 
sense  that  despinning  and  detumbling  can  be  accomplished  with  relatively 
small  control  requirements,  and  that  no  violent  dynamical  behavior  occurs 
even  when  internal  control  torques  only  are  applied.  Also,  the  use  of 
momentum  wheels  to  reduce  thruster  requirements  is  possible,  although  no 
evidence  is  presented  that  the  wheels  provide  any  major  benefit  in  terms 
of  efficient  system  control. 
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Fig.  3.  OMV  Thrust  Torques,  Momentum  Wheels  Uncoupled. 
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Fig.  4.  Gimbal  (Universal  Joint)  Torques,  Uncoupled  Momentum  Wheels 
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Fig.  8.  Angular  Velocities  of  OMV  with  u,  and  uQ  Feedback  Only 
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Fig.  9.  Angular  Velocity  <*>OJof  OMV  with  u^  and  Ug  Feedback  Only 
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Fig.  12.  Target  Spin  and  Precession  Angle  Rates,  u7  and  u„ 
Feedback  Only 


,  rad/sec 
<oot ,  rad /sec 


0.0  50.0  100.0  150.0  200.0  250.0  300.0 

TIME,  SEC 


Fig.  14.  OMV  Angular  Velocity  Components;  Full  Feedback  Added 
at  t  =  250  seconds 
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Fig.  15.  OMV  b3  Angular  Velocity  Component;  Full  Feedback  Added 
at  t  -  250  seconds 
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Fig.  17.  Gimbal  Control  Torque  u7  With  Full  Feedback  Added 
at  t  =  250  seconds  ' 
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Fig.  20.  OMV  Thruster  Torques,  b-  Momentum  Wheel  Torque 
Coupled  to  Target  Precession  Angle 
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Fig.  21.  t>2  Momentum  Wheel  Control  Torque,  u,-,  and  Wheel  Angular 

c  Vel.;  Control  Torque  u.  Coupled  to  Precession  Angle 
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TABLE 

III. 

G 

time 

U7 

U8 

TC 

|FC| 

FC 

2 

0 

•1.500 

-.192 

.000 

.613 

-.576 

5 

-.822 

-.242 

-.984 

.503 

-.225 

*• 

10 

-.176 

-.267 

-1 .576 

.631 

.097 

15 

.399 

-.274 

-1.897 

.759 

.374 

ft 

20 

.884 

-.267 

-2.032 

.859 

.601 

25 

1.273 

-.251 

-2.043 

.935 

.778 

30 

1.568 

-.230 

-1.973 

.993 

.907 

35 

1.777 

-.207 

-1.853 

1.033 

.992 

40 

1.909 

-.183 

-1.705 

1.056 

1.041 

45 

1.976 

-.161 

-1.545 

1.061 

1  .058 

50 

1.989 

-.140 

-1.382 

1.051 

1.050 

9 

55 

1.959 

-.121 

-1.224 

1.027 

1.023 

a 

60 

1.895 

-.104 

-1.076 

.990 

.981 

65 

1.808 

-.090 

-.939 

.945 

.928 

70 

1.703 

-.077 

-.815 

.892 

.869 

r. 

75 

1.588 

-.067 

-.704 

.835 

.806 

80 

1.468 

-.058 

— .  606 

.775 

.741 

85 

1.346 

-.050 

-.520 

.714 

.676 

90 

1.226 

-.043 

-.445 

.654 

.613 

* 

95 

1 . 109 

-.038 

-.380 

.595 

.553 

100 

.998 

-.033 

-.324 

.538 

.496 

105 

.894 

-.029 

-.276 

.484 

.443 

is 

110 

.797 

-.025 

-.235 

.434 

.393 

115 

.708 

-.022 

-.200 

.387 

.348 

120 

.626 

- .  020 

-.170 

.344 

.307 

•V 

125 

.552 

-.018 

-.144 

.304 

.270 

$ 

130 

.485 

-.016 

-.123 

.268 

.236 

135 

.425 

-.014 

-.104 

.236 

.207 

til 

140 

.371 

-.012 

-.089 

.207 

.  180 

1 

145 

.323 

-.011 

-.076 

.  181 

.156 

fjl 

150 

.281 

-.010 

-.064 

.157 

.136 

155 

.244 

-.009 

-.055 

.137 

.117 

h?' 

160 

.211 

-.008 

-.047 

.119 

.  101 

C* 

165 

.182 

-.007 

-.040 

.  103 

.087 

170 

.157 

-.007 

-.034 

.089 

.075 

175 

.135 

-.006 

-  .029 

.077 

.064 

s 

180 

.116 

-.005 

-.025 

.066 

.055 

185 

.100 

-.005 

-.022 

.057 

.047 

190 

.085 

-.004 

-.019 

.049 

.040 

5*3 

195 

.073 

-.004 

-.016 

.042 

.034 

S> 

200 

.062 

-.004 

-.014 

,035 

.029 

205 

.053 

-.003 

-.012 

.030 

.025 

- 

210 

.045 

-.003 

-.010 

.026 

.021 

_  _ 

215 

.038 

-.003 

-.009 

.022 

.018 

220 

.033 

-.002 

-.008 

.019 

.015 

225 

.028 

-.002 

- .  007 

.016 

.013 

t 

230 

.023 

-.002 

-  ,006 

.013 

.011 

•  *. 
s 

4  9 

#  % 
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TABLE  III., 

cont ' d 

time 

U7 

u8 

TC 

!FC| 

235 

.020 

-.002 

-  .005 

.011 

240 

.017 

-.002 

-  .004 

.010 

245 

.014 

-.001 

-.004 

.008 

250 

.012 

-.001 

-.003 

.007 

255 

.010 

-.001 

-.003 

.006 

260 

.008 

-.001 

-.002 

.005 

265 

.007 

-.001 

-.002 

.004 

270 

.006 

-.001 

-.002 

.003 

275 

.005 

-.001 

-.002 

.003 

280 

.004 

-.001 

-.001 

.002 

285 

.004 

-.001 

-.001 

.002 

290 

.003 

-.001 

-.001 

.002 

295 

.002 

-.001 

-.001 

.001 

7 

?? 


.00 

.00! 

.00* 

.00! 

.00! 

.00 

.oo; 

,oo; 

.oo: 

.oo: 

.oo: 

.oo: 

.oo: 


Note:  This  table  of  universal  joint  control  torques  and 
constraint  loads  applies  to  the  case  of  feedback 
control  with  uncoupled  momentum  wheels. 


TABLE  IV 


time 

U7 

U8 

TC 

1 FC  | 

Fc 

0 

-1.500 

-.192 

.000 

.613 

-.576 

5 

-.822 

-.242 

-.984 

.503 

-.225 

10 

-.176 

-.267 

-1.576 

.631 

.097 

15 

.399 

- .  274 

-1 .897 

.759 

.374 

20 

.884 

-.267 

-2.032 

.859 

.601 

25 

1.273 

-.251 

-2.043 

.935 

.778 

30 

1 .568 

-.230 

-1.973 

.993 

.907 

35 

1.777 

-.207 

-1.853 

1.033 

.992 

40 

1.909 

-.183 

-1.705 

1 . 056 

1.041 

45 

1.976 

-.161 

-1.545 

1.061 

1.058 

50 

1.989 

-.140 

-1.382 

1.051 

1.050 

55 

1.959 

-.121 

-1.224 

1.027 

1.023 

60 

1.895 

-.104 

-1.076 

.990 

.981 

65 

1.808 

-.090 

-.939 

.945 

.928 

70 

1.703 

-.077 

-.815 

.892 

.869 

75 

1.588 

-.067 

-.704 

.835 

.806 

80 

1.468 

-.058 

-.606 

.775 

.741 

85 

1.346 

-.050 

-.520 

.714 

.676 

90 

1.226 

-.043 

-.445 

.654 

.613 

95 

1.109 

-.038 

-.380 

,595 

.553 

100 

.998 

-.033 

-.324 

.538 

.496 

105 

.894 

-.029 

-.276 

.484 

.443 

110 

.797 

-.025 

-.235 

,434 

.393 

115 

.708 

-.022 

-.200 

.387 

.348 

120 

.626 

-.020 

-.170 

.344 

.307 

125 

.552 

-.018 

-.144 

.304 

.270 

130 

.485 

-.016 

-.123 

.268 

.236 

135 

.425 

-.014 

-.104 

.236 

.207 

140 

.371 

-.012 

-.089 

.207 

.  180 

145 

.323 

-.011 

-.076 

.181 

.156 

150 

.281 

-.010 

-.064 

.157 

,136 

155 

.244 

-.009 

-.055 

.137 

.117 

160 

.211 

-.008 

-.047 

.119 

.101 

165 

.182 

-.007 

-.040 

.  103 

.087 

170 

.157 

-.007 

-.034 

.089 

.075 

175 

.135 

-.006 

-.029 

.077 

.064 

180 

.116 

-.005 

-.025 

.066 

♦  055 

185 

.100 

-.005 

-.022 

.057 

.047 

190 

.085 

-.004 

-.019 

.049 

.040 

195 

.073 

-.004 

-.016 

.042 

.034 

200 

.062 

-.004 

-.014 

.035 

.029 

205 

.053 

-.003 

-.012 

.030 

.025 

210 

.045 

-.003 

-.010 

.026 

.021 

215 

.038 

-.003 

-.009 

.022 

.018 

220 

.033 

-.002 

-.008 

.019 

.015 

225 

.028 

-.002 

-.007 

.016 

.013 

230 

.023 

-.002 

-.006 

.013 

.011 

38 


TABLE  IV.  cont'd 


time 

U7 

U8 

TC 

1 FC  | 

Fc 

235 

.020 

-.002 

- .  005 

.011 

.009 

240 

.017 

-.002 

-.004 

.010 

.008 

245 

.014 

-.001 

-.004 

.008 

.006 

250 

.012 

-.001 

-.003 

.007 

.005 

255 

.010 

-.001 

-.003 

.006 

.005 

260 

.008 

-.001 

-.002 

.005 

.004 

265 

.007 

-.001 

-.002 

.004 

.003 

270 

.006 

-.001 

-.002 

.003 

.003 

275 

.005 

-.001 

-.002 

.003 

.002 

280 

.004 

-.001 

-.001 

.002 

.002 

285 

.004 

-.001 

-.001 

.002 

.002 

290 

.003 

-.001 

-.001 

.002 

.001 

295 

7? 

.002 

-.001 

-.001 

.001 

.001 

Note:  This  table  of  universal  joint  control  torques  and 
constraint  loads  applies  to  the  case  of  feedback 
control  with  the  b?  momentum  wheel  coupled. 


TABLE  V 


time 

U1 

U2 

U3 

U5 

2 

0 

-1.732 

-.066 

.  180 

.000 

.000 

5 

-.960 

-1.209 

.182 

.096 

.005 

10 

-.242 

-1.906 

.205 

.157 

.016 

15 

.382 

-2.296 

.238 

.  191 

.032 

20 

.893 

-2.473 

.269 

.206 

.051 

25 

1.286 

-2.504 

.291 

.207 

.069 

30 

1.566 

-2.438 

.301 

.199 

.088 

35 

1.744 

-2.311 

.299 

.  184 

.  105 

40 

1.832 

-2.147 

.287 

.  165 

.  121 

45 

1  .845 

-1.964 

.268 

.  145 

.135 

50 

1.799 

-1.777 

.244 

.123 

.148 

55 

1.706 

-1.593 

.218 

.102 

.158 

60 

1.579 

-1.418 

.192 

.082 

.166 

65 

1.430 

-1.256 

.167 

.062 

.173 

70 

1.267 

-1.108 

.144 

.045 

.178 

75 

1  .098 

-.974 

.122 

.028 

.181 

80 

.929 

-.855 

.104 

.014 

.183 

85 

,  .764 

-.750 

.088 

.001 

.183 

90 

,607 

-.657 

.074 

-.010 

.183 

95 

.461 

-.576 

.062 

-.019 

.182 

100 

.325 

-.506 

.052 

-.028 

.  180 

105 

.203 

-.445 

.044 

-.035 

.177 

110 

.093 

-.392 

.037 

-.040 

.173 

115 

-.005 

-.346 

.032 

-.045 

.  169 

120 

-.090 

-.306 

.027 

-.048 

.165 

125 

-.163 

-.271 

.023 

-.051 

.161 

130 

-.226 

-.242 

.020 

-.053 

.  156 

135 

-.278 

-.216 

.017 

-.055 

.151 

140 

-.322 

-.194 

.015 

-.056 

.146 

145 

-.357 

-.174 

.013 

-.056 

.141 

150 

-.385 

-.157 

.011 

-.056 

.  136 

155 

-.406 

-.142 

.010 

-.056 

,131 

160 

-.422 

-.129 

.009 

-.055 

.126 

165 

-.433 

-.118 

.008 

-.054 

.121 

170 

-.439 

-.108 

♦  007 

-.053 

.116 

175 

-.441 

-.099 

.006 

-.052 

.111 

180 

-.441 

-.091 

.006 

-.051 

.106 

185 

-.437 

-.084 

.005 

-.049 

.  102 

190 

-.432 

-.077 

.004 

-.048 

.097 

195 

-.424 

-.072 

.004 

-.046 

.093 

200 

-.416 

-.067 

.004 

-.045 

.089 

205 

-.405 

-.062 

.003 

-.043 

.085 

210 

-.394 

-.058 

.003 

-.041 

.081 

215 

-.383 

-.054 

.003 

-.040 

.077 

220 

-.371 

-.050 

.002 

-.038 

.074 

225 

-.358 

-.047 

.002 

-.037 

.070 

230 

-.345 

-.044 

.002 

-.035 

.067 

TABLE  V.  cont'd 


time 

U1 

U2 

U3 

U5 

2 

235 

-.333 

-.041 

.002 

-.034 

.064 

240 

-.320 

-.039 

.002 

-.032 

.061 

245 

-.307 

-.037 

.001 

-.031 

.058 

250 

-.295 

-.035 

.001 

-.030 

.055 

255 

-.283 

-.033 

.001 

-.028 

.053 

260 

-.271 

-.031 

.001 

-.027 

.050 

265 

-.259 

-.029 

.001 

-.026 

.048 

270 

-.248 

-.027 

.001 

-.025 

.046 

275 

-.237 

-.026 

.001 

-.023 

.043 

280 

-.226 

-.024 

.001 

-.022 

.041 

285 

-.216 

-.023 

.001 

-.021 

.039 

290 

-.206 

-.022 

.001 

-.020 

.037 

295 

T? 

-.197 

-.021 

.001 

-  .019 

.036 

Note:  This  table  of  thruster  torques  u,,  u2,  and  u3  and 
wheel  torque  ug  applies  to  the  cdse  6f  feedback 
control  with  5  the  b2  momentum  wheel  coupled. 
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TABLE  VI 


time  u7  u8  TC  |FC| 


■\) 

0 

-1 .500 

-.192 

.021 

.612 

-.500 

'0 

;  f  i 

5 

-.330 

-.281 

-.471 

.340 

-.225 

10 

-.403 

-.314 

-.738 

.242 

-.094 

tm 

15 

-.201 

-.298 

-.882 

.223 

-.071 

1 

20 

-.180 

-.251 

-.974 

.237 

-.124 

,v« 

25 

-.285 

-.192 

-1.054 

.288 

-.223 

30 

-.463 

-.136 

-1.139 

.379 

-.346 

>3 

35 

-.671 

-.093 

-1.231 

,495 

-.476 

& 

40 

-.878 

-.068 

-1.323 

.619 

-.602 

45 

-1.067 

-.063 

-1.402 

.740 

-.721 

50 

-1.230 

-.074 

-1.454 

.851 

-.831 

5 

55 

-1.371 

-.095 

-1.470 

.953 

-.936 

ft. 

60 

-1.499 

-.119 

-1.444 

1.051 

-1.038 

65 

-1.629 

-.139 

-1.379 

1.150 

-1 .141 

70 

-1.771 

-.151 

-1.280 

1.253 

-1.248 

►  *' 

VI. 

75 

-1.933 

-.151 

-1 . 159 

1.363 

-1.359 

80 

-2.119 

-.141 

-1.029 

1.478 

-1.473 

85 

-2.322 

-.122 

-.900 

1.596 

-1.589 

90 

-2.537 

-.099 

-.784 

1.712 

-1.704 

r 

95 

-2.751 

-.075 

-.687 

1.822 

-1.813 

100 

-2.954 

-.054 

-.612 

1.922 

-1.914 

/V.* 

105 

-3.139 

-.036 

-.559 

2.012 

-2.003 

c 

110 

-3.300 

- .  024 

-.524 

2.088 

-2.080 

115 

-3.433 

-.016 

-.504 

2.152 

-2.143 

KA 

M 

120 

-3.541 

-.011 

-.492 

2.204 

-2.194 

125 

-3.626 

-.010 

-.483 

2.245 

-2,235 

Hi 

130 

-3.693 

-.010 

-.473 

2.279 

-2,268 

135 

-3.746 

-.010 

-.458 

2.306 

-2,295 

i 

140 

-3.792 

-.010 

-.440 

2.330 

-2.318 

P 

145 

-3.834 

-.009 

-.416 

2,351 

-2,338 

150 

-3.873 

-.008 

-.391 

2.371 

-2.357 

155 

-3.912 

-.005 

-.364 

2.389 

-2,375 

*  M 

160 

-3.950 

-.003 

-.340 

2.407 

-2.392 

»V 

165 

-3.986 

.000 

-.318 

2.423 

-2.408 

170 

-4.020 

.003 

-.301 

2.437 

-2.422 

H? 

175 

-4.049 

.005 

- ,  288 

2.449 

-2.434 

;-:j 

180 

-4.072 

.006 

-.280 

2.458 

-2.443 

185 

-4.090 

.007 

-.276 

2,465 

-2.449 

< 

190 

-4.103 

.007 

-.274 

2.469 

-2.453 

,V 

■*ji 

195 

-4.110 

.006 

-.273 

2,471 

-2.455 

200 

-4.113 

.006 

-.272 

2.471 

-2.455 

205 

-4,114 

.005 

-.271 

2.470 

-2.454 

V 

#V 

210 

-4.112 

.004 

-.269 

2.469 

-2.452 

215 

-4,110 

.003 

-.266 

2,466 

-2.450 

220 

-4.107 

.003 

-.262 

2.464 

-2.447 

.  « 

225 

-4.104 

.003 

-.258 

2.461 

-2,445 

230 

-4.101 

.003 

-.253 

2,459 

-2.442 

L 
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TABLE  VI.  cont'd 


ime 

u7 

U8 

TC 

|FC| 

FC 

2 

35 

-4.099 

.003 

-.249 

2.456 

-2.440 

40 

-4.096 

.003 

-.246 

2.453 

-2.437 

45 

-4.093 

.003 

-.243 

2.450 

-2.434 

50 

-3.847 

.003 

-2.162 

2.495 

-2.284 

55 

-3.192 

.003 

-1.813 

2.072 

-1.896 

60 

-2.643 

.004 

-1.519 

1.720 

-1.572 

65 

-2.194 

.005 

-1.276 

1.432 

-1.306 

70 

-1.824 

.005 

-1.075 

1 .194 

-1.087 

75 

-1.519 

.006 

-.907 

.998 

-.906 

80 

-1.267 

.006 

-.768 

.836 

-.757 

85 

-1.059 

.006 

-.651 

.702 

-.634 

90 

-.887 

.006 

-.553 

.590 

-.532 

95 

-? 

-.743 

.006 

-.471 

.497 

-.446 

Note:  This  table  of  universal  joint  control  torques  and 
constraint  loads  applies  to  the  case  of  feedback 
with  control  torques  u,  and  ug  only  until  t  =  250 
seconds,  at  which  time'the  complete  control  vector 
u  is  fed  back. 


TABLE  VII. 

cont'd 

8 

time 

tj 

01 

W02 

"03 

y4 

>5 

235 

-.001 

.010 

.101 

.000 

-.001 

240 

-.001 

.010 

.  101 

.000 

-.001 

\> 

245 

-.001 

.010 

.  101 

.000 

-.001 

250 

-.001 

.009 

.101 

.000 

-.001 

m 

255 

.000 

.008 

.101 

.000 

-.001 

? 

260 

,000 

.006 

.101 

,000 

-.001 

265 

.000 

.005 

.101 

.000 

-.001 

270 

.000 

.004 

.101 

.000 

-.001 

275 

.000 

.003 

.101 

.000 

-.001 

£ 

280 

.000 

.002 

.  101 

.000 

.000 

285 

.000 

.002 

.101 

.000 

.000 

290 

.000 

.002 

.  101 

.000 

.000 

l 

295 

.000 

.001 

.101 

.000 

.000 

;<£»  NOTE:  This  table  of  OMV  angular  velocity  components, 

it  target  spin  rate  and  precession  angle  rate  of 

change  is  for  the  case  of  feedback  with  n-,  and 
ug  control  vector  components  only,  until  ' 

«jj  °  t  =  250  seconds  when  control  is  with  the 

^  complete  u  vector. 

s 


V.  Conclusion 


A  nonlinear  feedback  control  law  was  developed  and  used  to  despin 
and  detumble  an  axially  symmetric  target  satellite  originally  in  steady 
spin  and  precession.  The  control  law  derivation  is  based  upon  Liapunov 
stability  theory,  and  ensures  the  global  asymptotic  stability  of  the 
final  spin-stabilized  equilibrium  state.  The  results  indicate  that  the 
system  is  well  behaved,  in  the  sense  that  changes  in  both  the  system 
state  and  in  the  control  torques  are  smooth  throughout  the  maneuver.  The 
control  torque  magnitudes  are  relatively  small,  and  no  extreme  loading  of 
the  connecting  joint  between  OMV  and  target  satellites  occurred.  The 
system  could  be  driven  very  close  to  the  spin-stabilized  state  using  the 
joint  control  torques  alone.  However,  a  residual  target  precession  angle 
remained  at  the  end  of  the  300  second  maneuver,  as  did  non-zero  b-j  and  b2 
OMV  velocity  components.  This  is  due  to  the  fact  that  after  approximately 
200  seconds  of  feedback  control,  the  matrix.  A,  of  Eq  (26)  is  nearly  dia¬ 
gonal,  and  as  a  result  the  control  torques  Uy  and  ug  are  coupled  strongly 
only  to  the  target  precession  angle  rate  and  spin  rate,  and  the  precession 
angle  (states  Xy,  xg,  and  xg),  all  of  which  have  very  small  values.  The 
available  control  torques  at  the  universal  joint  are  thus  insufficient  for 
the  reduction  of  target  precession  angle  at  any  appreciable  rate.  Imple¬ 
mentation  of  full  vector  control  (using  all  eight  control  torques)  at 
t  =  250  seconds  successfully  drove  the  system  to  the  spin-stabilized 
equilibrium. 

Although  the  OMV  thrust  torque  magnitude  could  be  reduced  by 
coupling  a  momentum  wheel  torque  to  the  target  precession  angle,  no  at¬ 
tempt  was  made  to  accomplish  detumbling  with  momentum  wneel  and  joint 


'V 


I 


& 


-v 


* 


A 


»  ^ 


torques  alone.  Any  attempt  made  to  couple  the  bg  momentum  wheel  with 
system  states  or  Xg  resulted  in  an  increase  in  at  least  one  thrust 

torque  profile,  with  no  obvious  positive  influence  on  system  behavior. 
A  follow  on  effort  might  concentrate  on  the  development  of  a 


■s 

reliable  technique  for 

$ 

matrix  terms,  based  on 

p 

be  made  to  perform  the 

47 
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